A Minimalist Turbulent Boundary Layer Model 
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We introduce an elementary model of a turbulent boundary layer over a flat surface, given as a 
vertical random distribution of spanwise Lamb-Oseen vortex configurations placed over a non-slip 
boundary condition line. We are able to reproduce several important features of realistic flows, such 
as the viscous and logarithmic boundary sublayers, and the general behavior of the first statistical 
moments (turbulent intensity, skewness and flatness) of the streamwise velocity fluctuations. As an 
application, we advance some heuristic considerations on the boundary layer underlying kinematics 
that could be associated with the phenomenon of drag reduction by polymers, finding a suggestive 
support from its experimental signatures. 

PACS numbers: 47.27.-1, 47.27.nb, 47.27.De 



Turbulent boundary layers have been a central topic 
of interest in fluid dynamic research for long years [l|, Q ■ 
Nevertheless their obvious technological importance, a 
satisfactory description of the physical mechanisms which 
underlie the boundary velocity fluctuations remains elu- 
sive to date. As the result of intensive computational and 
experimental efforts carried out mainly along the last two 
decades, it is by now clear that the turbulent boundary 
layer is the stage for the production and the complex 
interaction of coherent structures [3|, a fact that was for- 
merly antecipated by Theodorsen '4| and Townsend 

Standard phenomenological formulations of the tur- 
bulent boundary layer problem aim at solving self- 
consistent equations for the expectation values of veloc- 
ity and the Reynolds stress tensor cornponents, relevant 
quantities in engineering applications [3, @| . At the very 
conception of these models no fundamental role is given 
to the whole boundary layer zoo of coherent structures, 
like streamwise and hairpin vortices, low speed streaks, 
etc. It is an open question, for instance, if a structural 
derivation of the Prandtl von-Karman logarithmic law 
of the wall is viable. In this sense, turbulent bound- 
ary layer modelling is a difficult problem of statistical 
physics, analogous to the derivation of thermodynamic 
equations of state from molecular kinematics/dynamics. 
The literature on the subject is still relatively small, al- 
though an initial discussion may be traced back to 1982 
with Perry and Chong [rj. 

This work differs from previous attempts 0, [^, [Tol . 
[TTI.IT^ essentially in its stronger simplifying assumptions, 
specific coherent structure modelling (the vortex-dipole 
model to be introduced below), and the analysis of higher 
order statistics for the streamwise velocity fluctuations. 
We do not seek, at the present level of mathematical 
treatment, quantitative agreement with experiments; in- 
stead, we compute a set of general profiles, which turn 
out to be clearly supported by observations. 

Our focus is on the streamwise fluctuations of the ve- 
locity field. Let us assume that close enough to the wall 
these fluctuations are mostly due to the flow generated 
by hairpin vortices [l3|, [HI, |l^ , like the one depicted in 
Fig.l, momentarily located in the surroundings of the 



measurement position. The main contribution to stream- 
wise fluctuations would come from the spanwise sector of 
hairpin vortices (also called "hairpin's head" ), while sub- 
dominant contributions would be related to their necks 
and legs. 
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FIG. f : A hairpin vortex which propagates along the posi- 
tive X direction. Streamwise fluctuations of the velocity field, 
associated to the passage of the hairpin vortex, are due es- 
sentially to the flow generated by its spanwise sector, which 
lies between points A and B. The curved arrow indicates the 
vorticity orientation. 



Since streamwise fluctuations of velocity are the only 
ones we have in mind, it is natural to replace the hairpin 
vortex by a simpler and more mathematically tractable 
structure, which we take to be a spanwise Lamb-Oseen 
vortex, an exact and non-stationary solution of the 
Navier-Stokes equations, described in cilindrical coordi- 
nates as 
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where F is the total circulation around the vortex and 
= 2vt is the squared vortex core radius at time t, 
defined in terms of the kinematical viscosity v. 
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Of course, ((T)) solves the fluid equations of motion in 
the absence of boundaries, so ^ is just a rough approxi- 
mation to a real vortex parallel to the wall. In our mod- 
elling definitions, we postulate that the symmetry axis 
of the vortex lies in the plane (henceforth designed the 
"measurement plane") that contains the measurement 
point and is normal to the wall. We assume, then, that 
at equally spaced time intervals, a given vortex is re- 
placed by another one, at a random distance y from the 
wall, with some prescribed probability distribution func- 
tion (pdf) p{y). 
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FIG. 2: The vortex dipole construction. The dashed line is 
contained in the measurement plane. The upper plane vor- 
tex is the "real" one, while the other is its mirror image. A 
uniform backgroung flow with velocity V is superimposed to 
the velocity field produced by the vortices, so that P be a 
stagnant point. 



As the model under consideration is effectively two- 
dimensional, a mirror vortex is introduced "below the 
wall" , so that streamlines do not cross the material sur- 
face. Furthermore, in order not to completely neglect 
the non-slip boundary condition, some improvement is 
attained if we make the velocity field to vanish at the 
intersection of the measurement plane with the wall. For 
this purpose, an external homogeneous velocity field is 
superimposed to the field generated by the vortex dipole. 
These definitions are shown in Fig. 2. 

An interpretation of the time dependence in Eq. ^ is 
in order. The time variable t is taken to be the total time 
elapsed since the hairpin vortex was created at the wall. 
This assumption, however, is not of great help, if there 
is no way to relate the vortex vertical position y to the 
time t. To solve this problem, at least in a phenomenolog- 
ical fashion, we find inspiration in the scaling structure 
of the laminar boundary layer over a flat surface. It is 
well-established, and analytically predicted by the Bla- 
sius solution, that the laminar boundary layer thickness 
grows with the distance from the leading edge as (5 ~ ^/x. 
This result can be understood in elementary terms as the 



fact that any small perturbation which is transported 
along the main direction of the flow (say, the horizontal 
one) follows a diffusive drift along the vertical direction. 
An analogy to the context of turbulent boundary lay- 
ers can be drawn, replacing words like "perturbations" 
by "coherent structures" and "molecular viscosity" by 
"eddy viscosity" . Actually, hairpin vortices are observed 
to grow in size as they get farther from the surface p^ . 
Therefore, we suppose that at time t, a diffusion-like re- 
lation t ^ y^ holds, and the vortex core radius can be 
written as rc = ay in where a plays the role of a 
phenomenological parameter in the turbulent boundary 
layer modelling 

The flow depicted in Fig. 2 is assumed to describe a 
boundary layer with no pressure gradient. This is so be- 
cause the velocity field is symmetric under reflections on 
the measurement plane. Pressure gradients are probably 
related to the flow induced by the whole gas of hairpin's 
vortices, appearing, therefore, as a "many-body" effect. 

We are now ready to work out a few relevant equations. 
Suppose that the Lamb-Oseen vortices are centered at y 
and —y. The streamwise velocity field vanishes at point 
P in Fig. 2. Once the resulting streamwise velocity is 
given as the sum of three contributions (the external, the 
real and the mirror vortex velocity fields), we get, using 
Eq. (P), 
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The above equation holds for any y only if the total cir- 
culation r is a y— dependent quantity. We find 
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It is a straightforward task to evaluate expectation values 
of general j/— dependent functionals F — F[u{y)] of the 
streamwise velocity field. It follows that 

{F[uiy)])= r dy'p{y')F[u{y,y')] , (4) 

where the "two-point velocity" , 



u(y,y') ^ V + 
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-[l-exp(-(y + y')72aV')] , (5) 



is nothing but the velocity field at y when the Lamb- 
Oseen vortices are placed at y' and —y'. 

We have applied (H]) for a set of velocity functionals, 
comparing the ?/— dependent profiles so obtained with ex- 
perimental results, as discussed below. As input param- 
eters, we take a = 1.0, V = 1.0 (i.e., the streamwise 
velocity is computed in units of the external velocity V). 
We use the pdf 



p{y) 



2b 



7r(y2 + 62) ' 



(6) 



3 



with b = 1.0, to model fluctuations of the vortex's height 
above the surface. It is important to note that the 
choice of the lorentzian distribution although arbi- 
trary, is by no means restrictive. We have checked out 
that smoothly decaying distributions lead to similar con- 
clusions, if one is indeed interested in a qualitative un- 
derstanding of turbulent boundary layer fluctuations. 

• The viscous and logarithmic layers 

The mean streamwise velocity is obtained from the ex- 
pectation value of F[u{y)] = u{y). The overall proflle is 
shown in Fig. 3, which interpolates between zero ve- 
locity at the wall and unit velocity at infinity. Even 
though u{y) seems to give a reasonable profile for the 
interval < y < oo, the model does not apply to the 
outer layer, because of the stronger interactions between 
coherent structures and also for the high intermittency 
produced by the random entrainmcnt of external laminar 
flow that take place in that region. 
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FIG. 4: The viscous layer, verified in the range 0.01 < y < 
0.25. The straight line has slope 0.95. 
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FIG. 3: The mean streamwise velocity. 



In Figs. 4 and 5, viscous and log-layers are clearly 
noticed for certain ranges of vertical distances. The ex- 
cellent fit of the data to the straight line in Fig. 5 is 
given by (u(y)) = 0.31 ln(j/) -I- 0.55. The numerical co- 
efficients have precision of 0.1%. Observe that a purely 
dimensional argument yields uiij) — Vf{y/b) in our par- 
ticular model. Therefore, the numerical verification of a 
log-layer forces us to identify the effective external veloc- 
ity V to the friction velocity, up to numerical factors. We 
may conjecture, thus, that the friction velocity is "what 
is left" when a few dominant vortical structures near the 
wall are removed. In other words, the friction velocity 
can be interpreted here as a "mean field" , while fluctua- 
tions are modeled by isolated vortex dipoles. 



FIG. 5: The logarithmic layer, verified for 0.6 < y < 1.85. 
The inset shows the same data plotted in linear scales. 



• Higher order statistics 

Let Fn[u{y)] = [u{y) — (u(y))]". We introduce 



Urms = V {F2[u{y)]) (7) 

and the hyperflatness functions Sn{y), given by 

{Fn[u{y)]) 



Sniy) 



(8) 



Fine measurements of turbulent boundary layer fluctu- 
ations for ^3 (skewness) and 54 (flatness or kurtosis), 
which can resolve the region very close to the surface, 
are reported in Ref. T?]. As we can see from Figs. 6, 7 
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and 8, there is a clear qualitative agreement with obser- 
vations. In particular, the abrupt sign-changing transi- 



tion of skewness is remarkably reproduced by the vortex 2.0- 
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FIG. 6: The Urms velocity. The inset shows the experimental 
measurements of by Lorkowski 




FIG. 7: The skewness for 0.3 < y < 5.0. The inset shows the 
experimental measurements of Ssd/)^^^ /urma by Lorkowski 
0- 



supressed at some point above the surface in turbulent 
boundary layers Also, by about the same time, in- 
teresting signatures of drag reduction have been found 
in the profiles of Urms{y) and S^{y) in connection with 
polymer dilution [2l[ . An additional peak is observed for 
Urms{y), while the skewness 83(1/) shows further sign- 
changing transitions. 




• Drag reduction by polymers 

The phenomenon of drag reduction by polymers (isl . 
[19I I is a long-standing puzzle of non-newtonian fluid me- 
chanics. The broad picture is that dissipation is attenu- 
ated due to the interaction of polymers with the coherent 
structures created near walls. 

Recent PIV (particle image velocimetry) experiments 
in flows with dilute polymers have shown that vortic- 
ity fluctuations - and probably coherent structures - are 



FIG. 9: The Urma velocity, as affected by coherent structure 
suppression, for 0.01 < y < 5.0. An additional peak is ob- 
served at y ~ 2.5. 



The vortex dipole model allows us to relate these ap- 
parently distinct features of drag reduction by polymers. 
Coherent structure supression can be naturally accounted 
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FIG. 10: The skewness, as affected by coherent structure sup- 
pression, for 0.2 < y < 2.1. Note the additional sign-changing 
transitions that take place within 1.5 < y < 2.0. 



for by changing the form of the vortex pdf p{y)- We 
take, for instance, a distribution which is uniform for 
< 2/ < c, but vanishes for y > c, that is. 



p{y)^c-'[eiy)-eiy-c)] , 



(9) 



where 0{y) is the Heaviside fmiction. Therefore, we sup- 
pose no vortex is found for y > c, as the result of polymer 
interactions. Setting c — 1.0 and keeping the previous 
definitions of a, 6 and V, we get the results shown in Figs. 



9 and 10, which are in striking correspondence with real 
profiles [2l|. 

We conclude with some general remarks. The present 
model does not take into account further aspects of the 
turbulent boundary layer phenomenology, which become 
important if the interest is shifted toward quantitative 
comparisons with experiments. An essential improve- 
ment, along the above lines, would be to introduce a pair 
of streamwise vortex configurations with opposite vortic- 
ity as a way to mimick hairpin legs. Only in this way it 
would be possible to compute the shear stress and, thus, 
define the physical scales of length and velocity which are 
necessary for the description of the inner boundary layer. 

To summarize, we have introduced an elementary 
model of a turbulent boundary layer, focusing our atten- 
tion on the streamwise fluctuations of velocity induced 
by hairpin vortices. The model's main scope is to pro- 
vide qualitative insights on the velocity and hyperflat- 
ness profiles. Up to the knowledge of the author, this 
is the first time the profiles of skewness and flatness of 
usual turbulent boundary layers, as well as certain statis- 
tical signatures of the phenomenon of drag reduction by 
polymers have been theoretically reproduced by means 
of vortex methods. 
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